In [1]:
library("dplyr")
library("readr")
In [2]:
medians <- readr::read_csv("NSQD_Res_TSS.csv") %>%
dplyr::group_by(parameter, units, season) %>%
summarize(median(res))
medians
Out[2]: